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Abstract 

Doppler laser cooling of a three-level ladder system using two near-resonant laser fields is analyzed 
in the case of the intermediate level being metastable while the upper level is short-lived. Analytical 
as well as numerical results for e.g. obtainable scattering rates and achievable temperatures are 
presented. When appropriate, comparisons with two-level single photon Doppler laser cooling 
is made. These results are relevant to recent experimental Doppler laser cooling investigations 
addressing intercombination lines in alkali-earth metal atoms and quadrupole transitions in alkali- 
earth metal ions. 
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I. INTRODUCTION 



Today, laser cooling of atoms is widely used in experiments where high precision spec- 
troscopy or precise control of the atomic motion is required. A large variety of schemes has 
been developed and applied to suit atoms with specific level structures and for obtaining 
particular temperature ranges In general, the lower the wanted temperature, the more 
sensitive does the light scattering process has to be on the velocity of the atom. For the most 
simple laser cooling scheme relying on the Doppler shift of an optical transition, Doppler 
laser cooling, this means the narrower the linewidth of the optical transition, the lower the 
obtainable temperature. However, since the maximum cooling force in the Doppler cooling 
scheme is dependent on the photon scattering rate, narrow linewidth transitions will lead to 
longer cooling times than wider transitions. 

The aim of the present study is to show that by driving, with two step-wise near-resonant 
laser fields, a three-level ladder system having an intermediate metastable level and a dipole 
emitting upper level, it is indeed possible to obtain both short cooling times and at the same 
time low temperatures through changes of the laser parameters. The present results concern 
in particular alkali-like ions Ba + , Sr + , or Ca + {2], but also contain information relevant for 
Doppler cooling of alkali-earth metal atoms using intercombination lines 3|, |4| . In addition 
to the excellent laser cooling properties, the presented scheme has the major advantage that 
individual ions can be detected by dipole emitted light from the upper level to the ground 
level at a wavelength being the sum of those of the two excitation fields. By the appropriate 
choice of laser parameters, a large scattering rate can be obtained, which can lead to an 
essentially background-free detection as it has been demonstrated experimentally recently 



Previously, ladder- (or cascade-) schemes have been studied in 0, 0, E, [8|] in the prospect 
of two-photon two-color Doppler cooling on narrow linos like demonstrated experimentally 
in jf| on alkaline-earth metals using an intercombination line. In these studies, the excited 
state of a laser-driven transition is coupled by a second laser to an upper more stable 
level. In the present study, the excited metastable level is coupled to a higher short-lived 
level, to enhance its coupling with the ground state. Such an excitation scheme is used 
for example for sideband cooling, but while sideband cooling processes have already been 
studied theoretically and experimentally 10, II ] in the low saturation regime, the purpose 
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of the present study is to go beyond the low saturation regime and to consider Doppler 
cooling in the regime where the sidebands are not resolved, as considered experimentally 
previously Q|. Our analysis applies to any atomic system that can be described by the 
ladder scheme depicted on Fig. [TJ Nevertheless, when numerical applications are needed, we 
refer to parameters describing the alkali-like ions Ba + , Sr + , or Ca + . 

This article is organized as follows. Sec. [Til introduces the theoretical model used to 
predict the population of the internal states. In Sec. II III the response of an atom at rest 
to the lasers excitation is studied, focussing on the expected signal. The Doppler cooling 
efficiency of the proposed excitation scheme is studied in Sec. [TV] and in Sec. [V] conclusions 
are reported. 



II. THEORETICAL MODEL 



The considered atomic system is composed by three electronic levels which are coupled 
by laser fields, according to the ladder-scheme depicted in Fig. [TJ State \g) is the ground 
state, \m) is metastable, and the excited state |e) decays radiatively into \g) and \m). 
This configuration of levels occurs, for instance, in alkaline-earth atoms such as Sr and Ca, 
rare-earth atoms like Yb and in alkali-like ions with a metastable <i-orbital, such as Hg + , 
Yb + , Ba + , Sr + , or Ca + . In these last cases, the relevant levels can be identified with the 
states \g) = 15*1/2), \m) = I-D3/2) and |e) = \P\/2)- Then, the transition \g) — > \m) is an 
electric quadrupole transition with a linewidth in the range from 0.01 to 0.1 Hz, while \g) 
and |m) couple to the excited state |e) with electric dipole transitions. We denote by 

= g,m,e) the resonance frequencies of the transitions. Transitions \g) — > |m), marked 
by w (weak), and |m) — > \e), marked by st (strong) are coupled by lasers at frequency uj w 
and u st , the strengths of the couplings are characterized by the Rabi frequencies Q w and Q st 
and the detunings defined as A w = uj w — u mg , and A st = u st — oj em . The radiative process 
coupling |e) to states \g) and \m) is characterized by the decay rate 7 and the branching 
ratio Peg/Pem (Pem + Peg = !)• The radiative decay of state \m), whose natural lifetime is 
of the order of 1 s for the mentioned ions, is neglected in the analytical model we present 
below. 

We denote by p the density matrix for the atom's internal degrees of freedom, while we 
treat the center of mass variables classically. x(t) denotes the time-dependent position of 
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FIG. 1: Ladder level scheme. The straight arrows represent the laser coupling between states \g), 
\m) and |e). The wavy lines indicate radiative decay taken into account in the model. 



the atom. The master equation for the density matrix p is 







IP 



[H,p]+Cp 



dt h' 

where Hamiltonian H gives the coherent dynamics and is decomposed into the terms 



H = Ho + H T 



where 



-hA w \m)(m\ - h(A w + A 8t )|e){e| 



gives the internal energies in the reference frames of the lasers while 

H T = ^e ifc ^|m)(5l + ^ £i e ifc ^«|e)(m|+H.c. 



(1) 



(2) 



(3) 



(4) 



2 1 /WI 2 

gives the laser coupling, where kj denotes the laser wave vectors (j = w, st). The relaxation 
operator 



= -7p(.p\ e )( e \ + |e)(e|p) 

+Pegl\g)(e\p\e)(g\ + (3 em -y\m){e\p\e)(m\ 



(5) 



describes the radiative processes. 

The response of the atomic system to laser excitation is revealed by the occupation prob- 
abilities of the atomic levels, which can be monitored by the number of scattered photons. 
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Detection of the atoms is made on the transition with the highest fluorescence rate |e) — > \g) 
(branching ratio j3 eg /(3 em > 1 for the ions mentioned above) which lies in the blue domain 
(< 500 nm). For the excitation scheme discussed in the following, this allows for photon 
detection without scattered light background since the excitation wavelengths lie in the red 
and near infra-red domain (> 680 nm). This major advantage does not hold for Hg + nor 
Yb + for which the excitation and emission wavelengths are very close. In the following, we 
focus on the case of Ba + , Sr + and Ca + for the choice of the relevant parameters needed to 
describe the atomic system. 

It should be noted that the dipole (strong) and the quadrupole (weak) couplings differ 
by some orders of magnitude. Nevertheless, powerful lasers in the red and near infra-red 
domain are now widely available and, assuming proper focussing of the laser beam, Rabi 
frequencies of the order of 1 MHz can be reached on the quadrupole transition. To give an 
idea of the feasible Rabi frequencies, let us mention that, for Ca + , 200 mW of laser power 
at 729 nm focussed on a waist radius of 50 fim results in Q w ~ 1 MHz [^j. 

III. INTERNAL DYNAMICS OF A LOCALIZED PARTICLE 

In this section we look for the solutions of Eq. ([1]) when the motion of the particle can 
be neglected. For numerical applications and when approximations are needed in analytical 
expressions, parameters leading to a strong driving of the |m) — > |e) transition with reduced 
sensitivity to Doppler shift are chosen { r y 2 <C Sl st , A 2 st ) . Furthermore, to remain close to 
experimental conditions, we do not require excessively large Rabi frequency on the \g) — > \m) 
transition and assume -C 7 2 , £l 2 st - 

In order to get some insight in the problem, the dressed state picture of the system 
can be of help in the limiting case where the weak coupling beam can be considered as a 
probe of the atom dressed by the photons of the strong transition. In this picture, the |m) 
state is broadened and shifted by the strong coupling Q st . For negative detuning A st we 
expect the \g) — > \m) transition resonance to be shifted to 5ls = — (y^st + + A st )/2. 
This is confirmed on Fig. [2] where the lineshape of the excitation probability in the |m) 
state p m = (m\p\m) shows a maximum very close to this value. Increasing the detuning 
A st decreases the light-shift as well as the power-broadened linewidth of the \g) — > |m) 
transition, resulting in an increase of the maximum of p m for constant coupling strength 
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FIG. 2: Steady state population of level |e) (a) and \m) (b) versus for 7/27r = 20 MHz, 
O st /2vr = 100 MHz and n w /2ir = 1 MHz. A st /2vr = -100 MHz (solid line), A S f/2vr = -200 MHz 
(dot-dashed line), A s t/2ir = —400 MHz (dashed line). In each case, the thick lines give result for 
Peg = 1 and the narrow line for f3 eg = 0.75. The vertical dashed lines show the light-shifts 5ls 
calculated for the three sets of curves. 

In an experiment, the measured signal is proportional to the probability p e = (e\p\e) to 
be in the |e) state. Its dependance with the detuning A st does not follow the same rule as 
p m and Fig. [2] clearly shows that a maximum excitation probability can be reached for a 
good choice of this detuning. This behaviour is studied in the following through resolution 
of the optical Bloch equations resulting from the master equation Eq. (CQ) 

A. Maximum fluorescence rate 

The steady state solution of Eq. (CQ) can be written analytically, assuming conservation 
of the total population (p g + p m + p e = 1) and without any approximation made. For the 
population p e , this solution can be written like 

Pe = n 2 st n 2 jv (6) 




6 



with 



V = (3 eg UA W (A W + A st ) UA W {A W + A st ) + 2Q 2 W - 2Q 2 st ) + 4A 2 y7 2 + Q 



+((3 eg - 1)20* (2A£ + 4A w A st - tf st ) + 2fi 2 (4(A U , + A si ) 2 + 7 2 + fi 2 ) . (7) 

This equation is used to estimate the maximum scattering rate that can be expected from 
this excitation scheme, the effective linewidth of the forbidden transition once coupled by 
the strong laser excitation, and the dependence of these two entities on the laser coupling 
parameters. Despite its complexity, Eq. [6] can be connected to the dressed state picture if 



one keeps in mind that for A^ = 5ls, 4A U) (A 1U + A st ) = Q 2 t . 

In the limit where Q 2 ^ <C j 2 \A w / A st \, the condition for maximum occupation probability 
of the |e) state, and thus for maximum fluorescence signal, can also be derived analytically 
from Eq. [6j In first order in 7 2 / (Q 2 t + A 2 J, this maximum is found for 

,2 \ 



mit + A 2 st) 



A w = 5 LS [1 - „ n2 ' ) ■ (8) 



This correction to the light-shift can be omitted in the evaluation of the maximum population 
p^ 1 as it adds a correction in second order in 7 2 /(f2 2 t + A 2 t ), which is negligible for the set 
of parameters of interest. For fi 2 <C Q 2 t , this maximum is 

p " = V I 4 + &7 + *• TSW + (& » " 11 S| + ?ft j ' (9) 

This last equation rules the maximum scattering rate, observed when the detuning of the 
weak transition is adjusted to the light-shift induced by the strong coupling. Assuming this 
detuning A w is adjusted, for each Rabi frequency Q w (which is experimentally the limiting 
factor) there is an optimal set (fl s t, A st ) which maximises p^ 1 . This optimal set can be found 
numerically but better understanding can be gained in the particular case of large detunings 
A 2 4 3> Q 2 t for which the light-shift can be approximated by the simple form 5ls — ^ 2 t/4A st . 
With this simplification, the optimal set is described by the analytic expression: 



which corresponds to the maximum population 

1 



l p ] f\ovt — / • (11) 

4 + 2(f3 eg - 1) + 2pfk gl ln w + 2 7 2 /fi 2 t 
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This expression is simplified in 



Pe I opt 



a. 



(12) 



7 2^2(3 eg 

if 7 2 /fi^ t is negligible compared to j/fiw 

Equation ( fl~2l) confirms the dependance of the maximum population p^' 1 with (3 eg , ob- 
served on Fig. [2J Indeed, this population increases when the branching factor /3 eg decreases. 
For the ions considered in this work, this factor is of the order of 0.75 for Ba + and 0.95 
for Ca + and Sr + . Since the previous equations show that different values of (3 eg result 
only in a variation of the maximum population (Eq. [T2l) and optimal set of parameters to 
reach it (Eq. [TO]) . (3 eg is fixed to 1 in the following. This simplification results mainly in a 
slight underestimation of the number of scattered photons but allows to handle more general 
equations. 
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FIG. 3: Steady state population of level |e) (a), |m) (b), \g) (c) versus — A s t for an atomic system 
close to Ba + , Sr + , or Ca + : 'y/2ir = 20 MHz, in the simplified case where (5 eg = 1. The different 
curves are for different values of Q w /2ir: 0.1 MHz (solid line), 0.5 MHz (dashed line) and 1 MHz 
(dotted line). O s t/27r = 100 MHz for all curves and A w is adjusted to 5ls for each value of A s t. 
The vertical dashed lines show the optimum detunings like defined by Eq. [TU1 for the three sets of 



curves. 
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Beside the dependance on the branching factors, it is worth noting that the optimal 
population given in Eq. [12] behaves like fi«,/7, and not like f^/7 2 , as a simple two level 
picture would have suggested. This linear behaviour is confirmed by numerical calculations, 
even out of the range of the simplifying approximation <C A^ t . This is illustrated on 
Fig. [3] where the probabilities of occupation of the three levels, derived from Eq. [T] with 
no approximations, are plotted for different values of Q w . For each detuning on the strong 
transition A S f, the weak transition detuning is adjusted to the light-shift due to the strong 
coupling: A w = 5l$- For a given value of Q st , these curves confirm an optimum value for 
A st which maximises the resonance value p^f . The numerical results also confirm the linear 
dependance of the maximal population with Q w and obey the simplified analytic solution 
given in Eq. [TJ] to better than 3%. Following Eq. [IU] the optimum of A st is expected to 
behave like 1 / y/£l w . This is also confirmed by the numerical results which show that the 
simplified equation Eq. [10] reproduces numerical observation for optimal A st to better than 
7% for the value of fl w chosen for this figure. 

We can compare this result with the case where the strong coupling is maximum (A st = 
0). Then Sls — ±^st/2 and in the regime described by fl s t > 7/2, the maximum of p^ 1 
is (Q w / , -f) 2 / j3 eg . For the realistic conditions we are concerned with (namely Q w < 7), the 
maximum population in |e) and thus the maximum number of scattered photons is smaller 
by a factor proportional to ^y/Q w to what can be expected for large detuning A st (see Eq. [12] 
for comparison). 

To know if the ladder-scheme can result in enough scattered photons to be efficiently 
detected, it is useful to compare this scheme with the one usually used for detection and 
Doppler cooling of the considered trapped ions 121 ] . It consists in the excitation on the 
\g) — > |e) and \m) — > |e) transitions, the last one being required for repumping ions de- 
cayed in the metastable state. This system forms a A-scheme and solving the optical Bloch 
equations for this system, one can show that the maximum occupation probability in the 
|e) state is close to 20 %. The parameters chosen for Fig. [3J show a maximum close to 2 
% for the ladder- scheme. The complete absence of background signal due to scattered light 
largely compensates a reduction of absolute fluorescence signal by an order of magnitude. 
Furthermore, an occupation probability of the excited state as high as 10 % can be found 
by numerical calculations with Rabi frequencies on the weak transition of 7 MHz. These are 
achieved by tight focussing of laser power superior to 500 mW. 
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B. Resonance linewidth 

The linewidth T e ff of the fluorescence resonance profile with the laser detuning A w and 
its dependance on the two couplings can be also studied numerically. This study gives 
evidence for two different types of behaviour for this linewidth (see Fig. HJ): When Q st is 
small compared to A st , the linewidth has a constant value, which increases linearly with 
Q w . Once Q st > A st , the linewidth grows with Q st and becomes independent on Q w . To 
reflect this behaviour in the analytic expression of the resonance linewidth, we first assume 
that Q w is the smallest parameter of the problem (f2^ <C fl 2 t , A 2 t , 7 2 , A^), to extract a limit 
linewidth we call T (FWHM): 



This expression is simplified in the regime where 7 2 <C Q 2 t + A 2 st and can be written like: 



Then, we assume that the increasing weak transition coupling power broadens this limit 
linewidth, resulting in the effective linewidth T e ff 



For strong saturation {VL 2 st / A 2 st ^> 1), Tq ~ 7/2(1 — A st /tt st ) tends to the maximum value of 
7/2, and T e ff ~ r . For low saturation, the limit linewidth behaves like r — 7/2(f2 2 t /2A 2 t ) 
and tends to 0, the main contribution to the linewidth is then due to the weak coupling: 

r e// ~ V2Q W . 

This dependance of T e jf with \Q st /A s t\ is shown on Fig. H]for different values of Q w and 
a far detuned strong coupling laser A st = — IO7. The comparison between Eq. [15] and the 
exact FWHM of the resonance profile of p e with A w , calculated numerically, shows a very 
good agreement as long as Q w remains weak compared to 7/2. When this approximation is 
not valid, the methods used above result in an overestimation of the linewidth for Q st > A st 
(dash-dotted line of Fig. HI). These curves show that the effective linewidth of the excitation 
scheme can be adjusted continuously by the ratio \Q s t/A s t\. When this ratio is chosen to 
optimize the maximum of p e (see Eq. [TUT) . the effective linewidth can be evaluated through 
Eq. [15j For VL W /^ <C 1, T — \/2Vt w and T e ff ~ 2VL W . This value gives the lower boundary 




(13) 




(14) 




(15) 
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FIG. 4: FWHM of the resonance curve of p e with A w , in units of 7/2, as a function of the ratio 
|r2 S f/A s t|, for a system identical to Fig. A s t/2ir = —200 MHz and f2 s j varies from 0.1|A s j| to 
10|A s j|. The different curves are for different values of £l w which are 0.1 MHz (solid line), 0.5 MHz 
(dashed line), 1 MHz (dotted line) and 5 MHz (dot-dashed line). The lines give the numerical 
results from the master equation, with no approximations. The symbols plots the FWHM as given 
by Eq. [El 



for an estimation of the effective linewidth when parameters are set to optimise the excited 
population. 



IV. TWO-STEP DOPPLER COOLING 



So far we have neglected the effect of the center-of-mass motion on the atomic dynamics. 
In the following, the considered atoms are trapped in a harmonic potential, like in Q , their 



motion is periodic and can give rise to sidebands [131 ) . To quantify the efficiency of the 
proposed two-step Doppler cooling, we assume that the effective linewidth V e ff is larger 
than the frequencies of oscillation Ui of the ions in the trap. Practically, it means that the 
sidebands are not resolved in the absorption spectra and 



rom the point of view of laser- 
atom interactions, the atoms can be considered as free [13|. Furthermore, we also assume 
;hat the recoil frequency hk^,/2m (x being w or st) is smaller than the effective linewidth 
14) 1 and that the frequencies of oscillation are high enough so that the root mean square 
of the velocity v rms is not significantly changed during one oscillation period. Within these 
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boundary conditions, the effect of the atom motion can be taken into account in the master 
equation (Eq. [I]) and its stationary state solution (Eq. [6]) by the Doppler effect which shifts 
the detuning from A st to A st — k st .v and A w to A w — k w .v. 

To simplify the notations, we consider a one dimensional problem, extrapolation to three 
dimensions being obvious. We assume that the two lasers propagate along the same axis, 
so v is the velocity along that axis and the positive sign is defined with respect to k w . 
Consequently, the projection k s t is positive (negative) if the laser on the strong transition 
copropagates (counterpropagates) with the laser on the weak transition. By definition, k w 
is always positive. 

For travelling waves, the radiative force F leading to Doppler cooling results from 
momentum transferred to the atom by photons on absorption/emission cycles. It is 
exactly the gradient of the atom-light interaction potential [14J] given in Eq. HJ F = 
hk w fl w lm(p gm ) + hk st Q st lm(p me ). Using the relations induced for the stationary state den- 
sity matrix by the master equation, one can show that this definition is equivalent to 
F(v) = Tik w f3 eg r yp e {v) + hkstjpeiv)- In the simplifying case where the decay from |e) to 
\m) is neglected {(3 eg = 1), the radiative force is 

F(v) = h(k w + k st ) We (v). (16) 

This is exactly the result one would obtain by analogy with a two level system where the 
momentum transferred on each absorption/emission cycle is h(k w + k st ) (when averaged on 
several cycles) and where the number of cycles per unit time is the excited state scattering 
rate 7 times the probability to be in this excited state p e . The introduction of the effective 
momentum \t e ff = (k^+k^) as the momentum transferred to the atoms is very significant. It 
demonstrates that the role of the coupling on the strong transition is not only to shorten the 
lifetime of the metastable state |m), which would be independent on the relative direction of 
the two lasers. This sensitivity of the cooling process to the relative direction of propagation 
has been experimentally demonstrated in 2J. 

F(v) can be split in a constant force F°, which acts on an immobile atom and in a velocity 
dependent force F(v) — F°, responsible for the cooling process. For trapped particles, F° 
shifts the center of the trajectory by F° /mujf and is not involved in the cooling. The cooling 
force, proportional to p e (v) — p e {v = 0) can be considered as a damping force — (3v only at 



first order in v {i.e. for v T e ff/k x ) [151]. 
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In the approximation where Q w is the smallest relevant coupling parameter, the analytic 
expression of the slope of p e with v is simple enough to give an insight to the dependence 
of j3 with the coupling parameters: 

k w + 8 (4A^(A st + A w ) — tt 2 st ) (A w (k st + k w ) + (A st + A w )k w 
(3 * -Hk w +k st ) 1Pe (0) Q%QUk w + k st ) 

(17) 

First, one can check that for negative detuning A st and negative relative detuning A w — 5ls, 
the j3 coefficient is positive and damping effectively takes place. Secondly, the particularity 
of the regime of large detunings A st and low saturation (f^ t <C A^ t ) for the strong transition 
is revealed by Eq. [T71 Indeed, in this limit, the detuning A w (which is very close to 6ls), is 
negligible compared to A st and then, only k w controls the slope of p e with v. It means that 
the velocity dependance of the force is proportional to k w v , and that, at first order in v, 
the efficiency of the strong coupling is not reduced by the Doppler effect. Nevertheless, the 
relative direction of the two lasers remains relevant for the effective momentum \a e ff trans- 
ferred to the atoms as the ratio between the sum and the difference of the two wavevectors 
k w and k st is equal to 11.9 for Ca + , 4.4 for Sr + and 1.9 for Ba + . 

When the light-shift becomes non-negligible compared to A st , the slope becomes steeper 
for copropagating lasers. However, the major effect of the relative direction of propagation 
remains the size of the effective momentum k e //. In the following, we consider only the 
case where the two lasers are copropagating and call k the maximum effective momentum. 
Fig. [5] shows the damping coefficient /3 for two different values of Q w , in a regime where 
the effective linewidth depends on Q w . The efficiency of the damping shows a maximum, 
which corresponds to a minimum of the cooling time, for detuning A„, — 5ls increasing with 
the linewidth of the signal. Numerical results show that, for fixed Rabi frequency on the 
weak transition Q w , (3 is maximised and so the cooling time minimised for Q st /A st close to 
its optimum defined to reach maximum fluorescence (see Eq. [TOl) . In practice, maximising 
the fluorescence rate also results in reducing the cooling time and this optimum condition 
was used for plots on Fig. [H The curves show that, with this condition, the maxima of the 
damping are very close and that the advantage of an increase in Q w is mainly to broaden 
the frequency range where damping is efficient. 

The temperature reached by Doppler cooling results from the equilibrium between the 
damping due to the momentum transfer and the diffusion induced by spontaneous emission. 
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FIG. 5: Damping factor /3 in units of hk 2 , versus the relative detuning of the weak transition from 
the light-shift. Ca + wavelengths are used (A,„ = 732 nm and X s t = 866 nm) and r y/2n = 20 MHz. 
The three curves are for different values of Q, w /2ir: 0.1 MHz (solid line), 0.5 MHz (dashed line) 
and 1 MHz (dotted line), /S. s t/2n is fixed to — IO7 and Sl s t is chosen to obtain the maximum signal, 
according to EqllOl 




By conservation of energy, the scattered photons have the momentum k and by analogy 



16) 



with two-level systems, the diffusion coefficient is 

D = h 2 k 2 1Pe (0) (18) 

and the temperature reached is ksT = D/ (3. 

Numerical calculations show that the minimum temperature Tp is reached for A w — 5ls — 
—T e ff/2 like expected by analogy with a two level system. They also show that the relevant 
parameter for the strong coupling is the ratio | O st / A st | since for identical ratio, different 
values of A st or Q st give the same results. The minimum temperature Tp is plotted on 
Fig. [6] versus \Q s t/A s t\, for the wavelengths of Ca + and with copropagating lasers. Like the 
effective linewidth, the limit temperature can be adjusted continuously down to few /xK by 
the choice of Q s t/A st , the lowest temperature being controlled by Q w . As a comparison, the 
limit temperature reached by Doppler cooling on the allowed dipole \g) — > |e) transition on 
Ca + ions is 0.55 mK. Contrary to a two-level system, the Doppler limit temperature can be 
lowered, but to the detriment of the fluorescence rate. Nevertheless, as low temperatures can 
be reached by continously reducing the effective linewidth, the cooling process can be started 
with a high damping factor, resulting in a total cooling time shorter than if the cooling 
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FIG. 6: Doppler limit temperature versus \Q s t/A s t\ for Ca + wavelength (X w = 732 nm and X s t = 
866 nm) and ^/2tt = 20 MHz, (3 eg = 1. For a given ratio \Q s t/A s t\ the curves are identical for 
different A s j. The different curves are for different values of Q w /2tt: 0.1 MHz (solid line), 0.5 MHz 
(dashed line) and 1 MHz (dotted line). 



process would be made permanently on a very narrow transition. As for the capture range, 
denned by the velocity for which the cooling force is maximal, numerical results confirm the 
intuitive idea that it increases with the effective linewidth T e ff, which in practice, increase 
with \Q 8 t/A 8 t\ and/or Q w (see part IIIIBj) . It is possible to reach a broad capture range at 
the same time as keeping a high damping factor and fluorescence rate by incresing \Q st /A st \ 
and obeying the optimum condition defined in lHI Al (Eq. ITOT) . This is shown on Fig. [7] where 
the dependence of the effective force F rms with the root mean square velocity is plotted 
for the condition leading to the minimum Doppler limit, A„, = 5ls ~ F e ff/2. Within the 
assumption defined at the beginning of this section, this effective force can be approximated 
by F(v rms ) - F(-v rms ). 

The lowest temperatures shown on Fig. [6] have only a theoretical meaning as, in practice, 
for reaching these temperatures, the effective linewidth of the system must be reduced 
and reach values smaller than the recoil and the oscillation frequencies. Then, to make 
predictions, other mechanisms than the Doppler cooling have to be taken into account, as 



the well known sideband cooling which wil 



has been experimentally demonstrated 



id, 



not be developed here but for which efficiency 
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FIG. 7: Effective force F(v rms ) in units of hk'y versus kv rms / (7 /2) for Ca~*~ wavelength (Au, — 732 
nm and X s t = 866 nm) and ^/2tt = 20 MHz, f3 eg = 1. The different curves are for different values 
of Q w /2tt: 0.1 MHz (solid line and dot-dashed line) and 1 MHz (dashed line). A s t is fixed to 
— IO7 and Q s t is chosen to obtain the maximum signal, according to Eq llOl only for the curves 
plotted with the solid line and the dashed line: Q s t = 34 MHz (solid line), 106 MHz (dashed and 
dot-dashed line). A w = 5ls — T e ff/2. 




V. CONCLUSIONS 



We have studied a two-step photoexcitation process in a ladder-system with a metastable 
intermediate state, as can be found e.g. in alkaline-earth atoms and in alkali-like ions. In 
case of Ca + , Sr + and Ba + , the two excitation wavelengths belong to the red or infrared 
domain whereas the scattered photons lie in the blue domain, which offers the possibility of 
fluorescence detection without stray light background. Our numerical and analytical studies 
show that this excitation process can be a powerful alternative for the detection of ions or 
readout of their internal state. Indeed, the number of photons can easily reached a tenth 
of what is expected from direct excitation on a dipole allowed transition and the absence of 
background signal due to scattered light can compensate this lower signal. Furthermore, we 
demonstrate that two-step Doppler cooling can efficiently take place and that resolution of 
oscillation sidebands is not required to reach temperatures lower than 100 fiK. The width 
of the resonance profile and the Doppler limit temperature can be varied continuously by 
tuning the power and/or detuning on the dipole allowed transition. This cooling mechanism 
can be of great interest for traps with low oscillation frequencies or with geometry leading 



1(3 



to a high stray light level like microtraps. 
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